Lie point symmetries of a general class of PDEs: The heat equation 
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■ Abstract 

We give two theorems which show that the Lie point and the Noether symmetries of a second-order 



ordinary differential equation of the form J ~ ^(^' (s) i (^)) ^'''^ subalgebras of the special 

projective and the homothetic algebra of the space respectively. We examine the possible extension of this 
result to partial differential equations (PDE) of the form A^-'uij — F(x^ ,u,Ui) = where u{x^) and Uij 
stands for the second partial derivative. We find that if the coefficients A^-' are independent of u{x^) then 
the Lie point symmetries of the PDE form a subgroup of the conformal symmetries of the metric defined 
by the coefficients A^-' . We specialize the study to linear forms of F{x^ ,u,Ui) and write the Lie symmetry 
conditions for this case. We apply this result to two cases. The wave equation in an inhomogeneous medium 



> 

oo 

fT^ ■ for which we derive the Lie symmetry vectors and check our results with those in the literature. Subsequently 



we consider the heat equation with a ffux in an n— dimensional Riemannian space and show that the Lie 
symmetry algebra is a subalgebra of the homothetic algebra of the space. We discuss this result in the case 
of de Sitter space time and in flat space. 
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X: 

- - ■ 1 Introduction 

A Lie point symmetry of an ordinary differential equation (ODE) is a point transformation in the space of 
variables which preserves the set of solutions of the ODE [1] [2l [3] . If we look at these solutions as curves in the 
space of variables, then we may equivalently consider a Lie point symmetry as a point transformation which 
preserves the set of solution curves. Applying this observation to the geodesic curves in a Riemannian (affine) 
space, we infer that the Lie point symmetries of the geodesic equations in any Riemannian (affine) space are 
the automorphisms which preserve the set of these curves. However, it is known by Differential Geometry that 
the point transformations of a Riemannian (affine) space which preserve the set of geodesies are the projective 
transformations. Therefore it is reasonable to expect a correspondence between the Lie symmetries of the 
geodesic equations and the projective algebra of the space. 
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The equation of geodesies in an arbitrary coordinate frame is a second-order ODE of the form 

F + rjfci^i'^ + F{x\x^) = (1) 

where F(x'' , x^) is an arbitrary function of its arguments and the functions r*^. are the connection coefficients of 
the spaced Equivalently equation ([l} is the equation of motion of a dynamical system moving in a Riemannian 
(affine) space under the action of a velocity dependent force. According to the above argument we expect that 
the Lie symmetries of the ODE ([Ij for a given function F{x'^, x^) are a subalgebra of the projective algebra of 
the space. This subalgebra is selected by means of certain constraint conditions which will involve geometric 
quantities of the space and the function F{x^ ,x^). This approach is not new and similar considerations can be 

found in gnnieiiaiiii]. 

The determination of the Lie point symmetries of a given system of ODEs consists of two steps (a) the 
determination of the conditions which the components of the Lie symmetry vectors must satisfy and (b) the 
solution of the system of these conditions. Step (a) is formal and it is outlined in e.g. [Tl[21|3]. The second step 
is the key one and, for example, in higher dimensions where one has a large number of simultaneous equations 
the solution can be quite involved. However, if we express the system of Lie symmetry conditions of (IT|) in a 
Riemannian space in terms of coUineation (i.e. symmetry) conditions of the metric, then the determination 
of Lie symmetries is transferred to the geometric problem of determining the generators of the projective 
group of the metric. In this field there is a vast amount of work that is already done to be used. Indeed 
the projective symmetries are already known in many cases or they can be determined by existing general 
theorems of Differential Geometry. For example the projective algebra and all its subalgebras are known for 
the spaces of constant curvature |10j and in particular for the flat spaces. This implies, for example, that the 
Lie symmetries of all Newtonian dynamical systems are "known" and the same applies to dynamical systems 
in Special Relativity! 

In this work we state a theorem which establishes the exact relation between the projective algebra of the 
space and the Lie symmetry algebra of ([1]), assuming that the function F depends only on the coordinates, i.e. 
F{x^). 

What has been said for the Lie point symmetries of ^ applies also to Noether symmetries (provided ([T]) 
follows from a Lagrangian). The Noether symmetries are Lie point symmetries which satisfy the constraint 

XmL + L^ = f. (2) 
dt dt ^ ' 

Noether symmetries form a closed subalgebra of the Lie symmetries algebra. In accordance to the above, this 

implies that the Noether symmetries will be related with a subalgebra of the projection algebra of the space 

where 'motion' occurs. As it will be shown, this subalgebra is contained in the homothetic algebra of the space. 

As it is well known, each Noether point symmetry is associated conserved current (i.e. first integral), hence 

the above imply that the (standard) conserved quantities of a dynamical system depend on the space it moves 

and the type of force F{x'^,x^) which modulates the motion. In particular, in 'free fall', i.e. in the case 

F{x^,x^) = the geometry of the space is the sole factor which determines the (standard) first integrals of 

motion. This conclusion is by no means trivial and shows the deep relation between Geometry and Physics! 



^This point of view can be generalized to the general second order ODE provided the functions F^j. can be identified with the 
connection coefficients of a metric. 
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A natural question which arises is the fohowing: 

To what extend this correspondence of Lie/ N aether symmetries of second order ODES of the form (QP, with 
the collineations of the space, is extendable to partial differential equations of second- order of a similar form? 

Obviously, a global answer to this question is not possible. However, it can be shown that for many 
interesting PDEs the Lie symmetries are indeed obtained from the collineations of the metric. Pioneering work 
in this direction is the work of Ibragimov [TT]. Recently, Bozhkov et al. [T^] studied the Lie and the Noether 
symmetries of the Poisson equation and shown that the Lie symmetries of the Poisson PDE are generated from 
the conformal algebra of the metric. 

In the present work we show that for a general class of PDEs of second-order, there is a close relation between 
the Lie symmetries and the conformal algebra of the space. Subsequently we apply these results to a number of 
interesting PDEs and regain existing results in a unified manner. As a new application, we determine the Lie 
symmetries of the heat equation with a flux in a rt— dimensional Riemannian space. 

The structure of the paper is as follows. In Section [2] we review briefly the collineations in a Riemannian 
space. In Section |3] we consider the equation of geodesies in an afhne space and determine the Lie symmetry 
conditions in covariant form. We find that the major symmetry condition relates the Lie symmetries with the 
special projective algebra of the space. A similar result has been obtained previously in [4] using the bundle 
formulation of second order ODEs. 

In Section 2] we solve, in a concise manner, the symmetry conditions and state Theorem [1] which gives the 
Lie symmetry vectors in terms of the collineations of the metric and Theorem [2] which gives the Noether point 
symmetries in terms of the homothetic algebra of the metric. 

In Section [5] we consider the PDEs of the form A^^Uij — F{x\u,Ui) = and derive the Lie symmetry 
conditions. We show that for these PDEs the Ci^'') provided A'^ /= and if A^(^ = then £^\x^)di is a 
Conformal Killing vector of the metric A^^ . 

In Section[6]we consider a linear form of F(x^, u, Ui) and determine the Lie symmetry conditions in geometric 
form. In section [7] we apply the results of section [B] to the wave equation in a two dimensional inhomogeneous 
medium and to the heat equation with flux in a n— dimensional Riemannian space. In the latter case, it is 
shown that the Lie symmetry vectors are obtained form the homothetic algebra of the A^^ metric. A similar 
result has been found for the Poisson equation [12]. 

Finally in Section |S] we comment on the results and point out possible directions for future research. 

2 Collineations of Riemannian spaces 

A coUineation in a Riemannian space is a vector field X which satisfies an equation of the form 

Cx^ = B (3) 

where Cx denotes Lie derivative, A is a geometric object (not necessarily a tensor) defined in terms of the 
metric and its derivatives (e.g. connection, Ricci tensor, curvature tensor etc.) and B is an arbitrary tensor 
with the same tensor indices as A. The collineations in a Riemannian space have been classified by Katzin et 
al. [13]. In the following we use only certain collineations. 
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A conformal Killing vector (CKV) is defined by the relation 

Cx9^J = {x^) 9^r (4) 

If ■(/; = 0, X is called a Killing vector (KV). If -0 is a nonvanishing constant then X is a homothetic vector (HV) 
and if ■i/' y = 0, X is a special conformal Killing vector (SCKV). A CKV is called proper if it is not a KV, a HV 
or a SCKV. 

A Projective collineation (PC) is defined by the equation 

CxT)^=2<j,^Jly (5) 

If (/) = 0, the PC is called an affine collineation (AC) and if ^.j^ = 0, a special projective collineation (SPC). A 
proper PC is a PC which is not an AC, HV or KV or SPC. The PCs form a Lie algebra whose ACs, HV and 
KVs form subalgebras. It has been shown that if a metric admits a SCKV, then also admits a SPC, a gradient 
HV and a gradient KV HiJ. 

In the following we shall need the symmetry algebra of spaces of constant curvature. In [TO] it has been 
shown that the PCs of a space of constant nonvanishing curvature consist of proper PCs and KVs only and if 
the space is flat then the algebra of the PCs consists of KVs/HV/ACs and SPCs. 

In particular, for the Euclidian space i?" the projection algebra consists of the vectors of in Table 1. 



Table 1: Collincations of Euclidean space i?" , I , J — 1,2, ... ,n. 



Collineation 


Gradient 


Nongradient 


Killing vectors (KV) 


S/ = S^d, 




nomothetic vector (HV) 


H = x'd. 




Affine CoUineation (AC) 


Ah = xiS]dt 


Aij = xj5]di 


Special Projective collineation (SPC) 




Pi = 5/H. 



3 The Lie point symmetry conditions in an affine space 

We consider the system of ODEs: 

n 

F + rjfci^i'^ + J2 Pji-j,^^' . . . i^" = (6) 

where F*,, are the connection coefficients of the space and P^_^ j^{t,x^) are smooth polynomials completely 
symmetric in the lower indices and derive the Lie point symmetry conditions in geometric form using the 
standard approach. Equation ([6]) is quite general and covers most of the standard cases, autonomous and non 
autonomous, and in particular equation ([T]). Furthermore because the rjj.'s are not assumed to be symmetric, 
the results are valid in a space with torsion. Obviously they hold in a Riemannian space provided that the 
connection coefficients are given in terms of the Christofell symbols. 

The detailed calculation has been given in [9] and shall not be repeated here. In the following we summarize 
these results. 
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The terms i^^ . . . i-'™ for m < 4 give the equations: 

L^P' + 2U + iP\t +V\u +V',t P) = (7) 

L^P] +ttP;+ ^P;,t + S] + 2^,, si) P*^ + 2rj\tu -U si + 27j\t P)k - (8) 

L,P;^ + L.rj-fe + S'l^ + t(k S\ai) P% + en,,t - 2e,t(, 51) + i^f^t P'dk, = (9) 

L^P'.kd - tt P)m + e,e Sl,P%) + CP^,,,, + 4,7^t P^,,, - ^(,^,,,5^) = (10) 

and the conditions due to the terms x^^ . . . i-'™ for m > 4 are given by the foUowing general formula: 
LvPl.,^+Pn..,^A + (:^-m)^,tPl.,^ + 

+ (2 - (m - 1)) P;,..,_/L + (- + 1) Pl..^..v!r' + e.^..,,„_/L = 0- (11) 
We note the appearance of the term Lri^^j^ in these expressions. 

Eqn ([T]) is obtained for m — Q, P' = P* in which case the Lie symmetry conditions read: 

L,,P' + 2UP' +^P\t+v\u=Q (12) 

S) + 2tj SI) P'= + 2r?\,|, -^,u SI = (13) 

L,T), - 2Uu SI) = (14) 

^U\kSd) = 0. (15) 

If P* = we obtain the Lie symmetry conditions for the geodesic equations (see [3]) ■ 

4 The autonomous dynamical system moving in a Riemannian space 

We 'solve' the Lie symmetry conditions (|12p - (IT5|) for an autonomous dynamical system in the sense that we 
express them in terms of the collineations of the metric. 

Equation (jisp means that ^ j is a gradient KV of gij. This implies that the metric gij is decomposable. 
Equation (jl4p means that rj^ is a projective coUineation of the metric with projective function ^ j. The remaining 
two equations are the constraint conditions, which relate the components ^, of the Lie symmetry vector with 
the vector F^{x^). Equation (IT2|) gives 

(P„5*^") P, + g'^L^F, + 2^^,g'iF, + = 0. (16) 

This equation is an additional restriction for rf because it relates it directly to the metric symmetries. Finally 
equation gives 

- 5)i,u + {S.,,51 + 2<5}C,fe) + 2,?:,^ + 2r},ry;; = 0. (17) 

We conclude that the Lie symmetry equations are equations . pT)) where ^(t, a;) is a gradient KV of the 
metric gij and rf {t, x) is a special Projective coUineation of the metric gij with projective function ^ j. We state 
the results in theorem [T] [5] . 
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Theorem 1 The Lie point symmetries of the system of equations of motion of an autonomous system under 
the action of the force (x* ) in a general Riemannian space with metric gij , namely 

x' + rjfei^i'^ = F' (18) 

are given in terms of the generators of the special projective algebra of the metric gij. 



If the force is derivable frora a potential V{x^) so that the equations of motion follow from the standard 
Lagrangian 

L(x^iJ) = ^g.jx'x^ -V{x^) (19) 



with Hamihonian 



the Noether conditions, are 



E = l-g^j£x^ + V (x^) (20) 



Vkv'^ + V^^t = -f,t (21) 

tt9^3-^,,y = (22) 

Ln9^j = 2 Q^,t^ 9i] (23) 

= 0. (24) 



Last equation impUes ^ = ^ (t) and reduces the system as fohows 



LriQij = 2 Q?,t^ 9i] (25) 
+ - -/,* (26) 

= (27) 

Equation ([25]) imphes that 77* is a conformal Kilhng vector of the metric provided t 7^ 0- Because gij is 
independent of t and ^ = £,{t) the rf must be is a HV of the metric. This means that rf {t,x) ^T{t)Y' (x^) 
where is a HV. If ^ ^ = then 77* is a KiUing vector of the metric. Equations (1^51) . ([Tf)) are the constraint 
conditions, which the Noether symmetry and the potential must satisfy for the former to be admitted. These 
lead to the fohowing theorem [9]. 



Theorem 2 The Noether point symmetries of the Lagrangian U9\) are generated from the homothetic algebra 
of the metric gij . 

More specifically, concerning the Noether symmetries, we have the following 

All autonomous systems admit the Noether symmetry dt whose Noether integral is the Hamiltonian E ([201) • 
For the rest of the Noether symmetries we consider the following cases 

Case I Noether point symmetries generated by the homothetic algebra. 
The Noether symmetry vector and the Noether function G {t, x*^) are 

X = 2'tljYtdt + Y'd^ , G {t, x^) = pt (28) 
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where, V'y is the homothetic factor of (V'y = for a KV and 1 for the HV) and p is a constant, provided 
the potential satisfies the condition 

CyV + 2ijYV +p^Q. (29) 

Case II Noether point symmetries generated by the gradient homothetic Lie algebra, i.e., both KVs and 
the HV are gradient. 

In this case the Noether symmetry vector and the Noether function are 

X = 2ijY Jt {t) dtdt + T {t) H'^^ , G {t, x^) = TtH {x'') + p J Tdt (30) 

where, is the gradient HV or a gradient KV, the function T{t) is computed from the relation T^tt — mT where 
TO is a constant and the potential satisfies the condition 

LhV + 2i;YV + mH +p = Q. (31) 

Concerning the Noether integrals we have the following result (not including the Hamiltonian) 

Corollary 3 The Noether integrals of Case I and Case II are respectively 

Icj = 2il;YtE - g.jY'x^ + pt (32) 

= 2^-^ jT{t)dt E- g^jH-'x^ + TtH + p J Tdt (33) 
where, E is the Hamiltonian i20\} . 

We remark that theorems [T] and [2] do not apply to generalized symmetries [151 [16]. 

In a number of recent papers [7j [171 (HI HS] , the authors study the relation between the Noether symmetries 
of the geodesic Lagrangian 

L = i.g.,i^i^' (34) 

where — {s is an affine parameter along the geodesies) with the spacetime symmetries. They also make a 
conjecture concerning the relation between the Noether symmetries and the conformal algebra of spacetime and 
concentrate especially on conformally flat spacetimes. In [TH] it is also claimed that the author has found new 
conserved quantities for spaces of different curvatures, which seem to be of nonnoetherian character. Obviously 
due to the above results (see also [8]) the conjecture/results in these papers should be revised and the word 
'conformal' should be replaced with the word 'homothetic'. 

It would be of interest to examine if the above close relation of the Lie and the Noether symmetries of the 
second order ODEs of the form ([6]) with the coUineations of the metric is possible to be carried over to some 
types of second order partial differential equations (PDEs) . Although to this question it is not possible to give 
a global answer, due to the complexity of the study and the great variety of PDEs, it is still possible to give 
an answer of some generality which concerns many interesting and important cases. We will do this in the 
remaining sections. 
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5 The case of the second-order PDE's 



In the attempt to establish a general relation between the Lie symmetries of a second order PDE and the 
coUineations of a Riemannian space we derive the Lie symmetry conditions for a second order PDE of the form 

A'^u,j - F{x\u,u,) = (35) 

and we consider the coefficients A^^ (x, u) to be the components of a metric in a Riemannian space. According 
to the standard approach [H El E] the symmetry condition is 

X^^\H)=\H (36) 

where A(a;*, u, Ui) is a function to be determined. X^^l is the second prolongation of the Lie symmetry vector 

^ = (37) 



given by the expression 

dx^ ' ''du ' du, ' '''J du. 



where 



(1) ^ ,3 .3 

(2) Drff' Z)^'" \ ffe , /^fc , tfe \ 

The introduction of the function X{x^ ,u,Ui) in p6p causes the variables x^^u^Ui to be independent!^. 
The symmetry condition X^'^\H) = \H when applied to ([35]) gives: 

A'^Tjf + {XA=) - (F) = X{A'u,, - ^^) (39) 

from which follows: 

= A'r^,^ - v^^g^'Fu, - X{F) + \F 

- A^ {UijUa + UiUja + UiaUj) - UiUjUaA^ Cuu- (40) 

We note that we cannot deduce the symmetry conditions before we select a specific form for the function F. 
However we may determine the conditions which are due to the second derivative of u because in these terms 
no F terms are involved. This observation significantly reduces the complexity of the remaining symmetry 
condition. Following this observation we have the condition: 

= A'w, - ^'^■(^';,u,a + tj^^k) + (e'4 + V^i^u,, - XA'^u,, 

- A^ {u^jUa + U^Uja + UiaUj) - UiUjUaA^ ^uu 



^See Ibragimov [11) p. 115 
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from which follow the equations: 

A''^ {uijUk + UjkU, + UikUj)^'^,^ = 

The first equation is written: 

^''^':,u + A'''C.,u + A'''e..u = ^ A^"'£!:l = o- (41) 

The second equation gives: 

A^'Vu + VA% + - A^^Ct - A^'^tk - A^'^' = 0. (42) 

and the last equation gives: 

A^'Cu = 0. (43) 
It is straightforward to sho-^J^I that condition (|4ip implies 



^tu = (44) 

which is a well known result. 

From the analysis so far we obtain the first result: 

Proposition 4 For all second-order PDEs of the form A^^Uij — F{x^ ,u,Ui) — 0, for which at least one of the 
A^^ is ^ Q the C u — ^ or = {x^ ) . Furthermore condition [JMl is identically satisfied. 



There remains the third symmetry condition (1421) . We consider the following cases: 

We write (|42| ) in an alternative form by considering A^^ to be a metric as follows: 

L^.o^A^' = XA^' - {7jA'^),u (45) 

from which follows: 

Proposition 5 For all second-order PDEs of the form A^^Uij — F{x^ ,u,Ui) = 0, for which A^^ ,u = i.e. 

A^^ = A'^{x^), the vector ^^di is a CKV of the metric A^^ with conformal factor (X — rj^){x). 



Assuming^ A** = A*' = we have 
for i = i = Q nothing 
for i, jj^O gives (|45|) and 
for i = 0, ^ becomes: 



A'^Vu + nA% + tAl - A'^^C* - A'\l^ ~ XA'^ = 

A'^Ck = (46) 



which leads to the following general result. 



^We give a simple proof for ra = 2 in Appendix A. A detailed and more general proof can be found in I20| . 
*The index t refers to the coordinate a;" whenever it is involved. 
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Proposition 6 For all second-order PDEs of the form A^^ Uij — F{x^ ,u,Ui) — for which A^^ is nondegenerate 
the ^* J, = 0, that is, ^* = . 

Using that = when at least one of the Aij ^ 0, the symmetry condition (pUj) is simplified as follows 

= A^^?7,, - ^^,A'F,u, - X{F) + AF 
+ lA'n^^u, ~ A'i%u^ - u^ii^A^F^^^ + e,^UkA'Fu, 

+ A*^'77„„^.,u, + A^r^^u,, ~ '^A^'£!:,^u,k (47) 
+ {^Al + 77^:^)^,, - XA^u,, 

which together with the condition P5|) are the complete set of conditions for all second-order PDEs of the 
form A-^Uij — F{x'^, u, Ui) — 0, for which at least one of the Aij ^ 0. This class of PDEs is quite general. This 
fact makes the above result very useful.. 

In order to continue we need to consider special forms for the function F{x^ u, ui). 



6 The Lie symmetry conditions for a linear function F{x^ u, ui) 

We consider the function F{x, u, Ui) to be linear in Wi, that is to be of the form 

F{x,u,Ui) ^ B''{x,u)uk + f{x,u) (48) 
where B^{x, u), f{x, u) are arbitrary functions of their arguments. In this case the PDE is of the form 

A'^u,, - B''{x, u)uk - fix, u) - 0. (49) 

The Lie symmetries of this type of PDEs have been studied previously by Ibragimov|ll). Assuming that at 
least one of the Aij ^ the Lie symmetry conditions are (|T7|) and (|^5|) . 
Replacing F{x,u,u\) in (|T7)) we fincij: 

= A^^,^ - n^.g'^B, - - r)f.u + A/ 
+ 1A^r]^,u, - A^i%na - u.r^^g'^B, + e';"fc.9*^B, + XB^^u^, - ,^B%Uk - ^'B'^^Uk 



+ A'^'ig^^u.uk (50) 

+ Ah^^u,, - 2A^^i\^u,, + (C'^4 + nA^^)u,, - XA'^u,, (51) 

from which follow the equations: 

A^i^,^ - 77,.S' - efM - Tlf,u + Xf = (52) 

-2A''Vm + + ^uB"" - e,.B' + Ci?" - AS'- + 77i?i = (53) 

A'^'Vuu = 0. (54) 



^We ignore the terms with Uij because we have aheady use them to obtain condition II45I I. Indeed it is easy to see that these 
terms give A^^rj^ — 2A^^£,'^- + C'^A'^ + rjA^^i — AA*-* = which is precisely condition l l45t . 
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Equation ([51]) gives (because at least one A^'' 7^ 0!) : 

T] = a{x')u + b{x'). (55) 

Equation (|53l) gives 

-2A'''a^., + aB'' + auB% + ^^^'^^^i^ - S^^-^B' + ^B':, - XB'' + bB% = 0. 

We summarize the above results as follows. 

The Lie symmetry conditions for the second order PDEs of the form 

A'^Uij - B'^{x, u)uk - f{x, m) = (56) 

where at least one of the Aij ^ are: 

A'^ {a,ju + h,) - {a^,u + b ,)B' - ^V^fc - auf^u - bf,,, + A/ = (57) 

A'^^% - 2A'''a^^ + aB^ + auB% - ^';.^B' + C B^, - \B^ + bB';^ = (58) 

L^.o^A^' = iX-a)A'' -vA'\u (59) 

ri = aix')u + b{x') (60) 
e'n-0^e'(^'). (61) 

We note that for all second order PDEs of the form A^^Uij — B^{x,u)uk — f{x,u) ~ for which A'-'.u = 
i.e. A'^^{x^), the C*(^"') is a CKV of the metric A'^. Also in this case A(a;'). This result establishes the relation 
between the Lie symmetries of this type of PDEs with the collineations of the metric defined by the coefhcients 
Aij. 

Furthermore in case the coordinates are t, x^ (where i = 1, n) A'* = A*^ = and A'^^ is a nondegenerate 
metric we have that 

e,^ = o^e{t)■ (62) 

These symmetry relations coincide with those given in |llj . 
Finally note that equation (|58|) can be written: 



A'3^% - 2A^^a,, + [C, Bf + (a - \)B^ + (aw + b)B)^ = 0. (63) 

Having derived the Lie symmetry conditions for the type of PDEs of the form A^^Uij ~ B^{x, u)uk — f{x, u) =0 
we continue with the computation of the Lie symmetries of some important PDEs of this form. 
Before we proceed, we state two Lemmas which will be used in the discussion of the examples. 



Lemma 7 a. In flat space (in which F^j, = 0) the following identity holds: 

i«r^-4- (64) 

b. For a general metric g-ij satisfying the condition L^ig.gij = — (A — a)gij the following relation holds: 

.]k -9 t., jk 



ff^%F:^., =5-'''C.fc+r''/-r,zr' + (a-A)r. (65) 
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Proof 

Using the formula 

I^i^^.jk — ^\jk,lt + jk ~ C ,l^^.jk + Cj^^.sk + C,k^]sj 

we have: 

= r:/+9"'e,,k-c,ir' 

- + .9^''e\ ,k - e\;r' + 2{g^'e,,T\i - 9''e,,^\i) - (A - a)P 

= rx'-r,;r' + ff^'^r,,fc + («-A)r 

Lemma 8 Assume that the vector is a CKV of the metric g^j with conformal factor —(A — a) i.e. L^ig.gij = 
— (A — a)gij. Then the following statement is true: 

5^%rV, = ^(a-A)-^ (66) 

where n = g-'^gkj ths dimension of the space. 
Proof 

Using the identity: 

L^^'jk = -jf [VfcL^gjr + VjL^gkr - VrL^gkj] (67) 
and replacing L^gij = (a — \)gij we find: 

Li^'jk = ^5"' [(a - A),fcgjr + (a - A)^^^^ - (a - A),^^^^] 
= ^ [(a - A),fc5} + (a - A) j^fc - g'^'ia ~ X)^r9kj] ■ 
Contracting with g^^ we obtain the reguired result. 

7 Applications 

7.1 The wave equation for an inhomogeneous medium 

In order to show how the above considerations are apphed in practice we consider the wave equation for an 
inhomogeneous medium in flat 2d Newtonian space: 

C^(xl)uil -M22 = 0. (68) 

In this case we have: 

An - c-2(xi), A22 = -1, A12 = 0; = 0; / = 0. 
The symmetry conditions (157)) - (pT|) become: 

A'^{aijU + h,j) = Q (69) 
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A'^i% - 2A*'=a. (70) 

L^,qA'' = {\-a)A'^ (71) 

T] = a{x')u + b{x^) (72) 

e{x). (73) 



The vector is a CKV of the metric Aij =diag(c~^(a;^), — 1) with conformal factor —(A — a). This is 
nondegenerate 2-d metric which is conformally flat, therefore if we find the conformal factor we will have the 
solution and the function a (a;*). 

We take now the metric to be the A^^- Then according to Lemma [7] we have (where F^^ j, are the V s of the 
metric Aij) : 

A^'^L^r^, = r^e - C/r' + A^'Cjk - (A - a)r ^ 

A^'^Cjk = A^'^Li^'.jk + - r:,e' + (A - a)r (74) 

Then the Lie symmetry condition ([7(I| A^^£_\j — 2A^^ai is written: 

A^^{L^T% - 26';a,) = T.^' - ^'^r' - (A - a)r\ (75) 

For the metric Ay — diag{~c'^{x^), —1) we compute F* = F'j = hence the last equation becomes: 

A^^(L^F^ -25^*a,) -0. (76) 

Because the metric Aij is nondegenerate this implies 

L^F^ = 2J^a,,) (77) 

which means that is a projective vector of the metric Aij with projection function a. But is also a CKV 
of the same metric with conformal factor a — A therefore must be a HV of the metric Aij . This implies 

a = ci a constant 

and 

-A + Ci = C2 ^ A = Ci - C2 

where C2 is the homothetic factor. From the remaining condition ((M)) we have A^^b,ij — that is, the function 
6 is a solution of the original wave equation We conclude that the Lie symmetry vector is (see also [T] p. 
182): 

X = C^^ + {ClU + b)du (78) 

where is a HV (not necessarily proper) of the metric Aij and 6(2:') is a solution of the wave equation. The 
vector S^^di is the sum UiKV^ + anvHV where ai,aHV are constants and KV-jHV are any of the KVs and 
the HV (if it exists) of the metric Aij . 

In the following section we consider a new example which is the heat conduction equation with a fLux in an 
n— dimensional Riemannian space. 



^Because det Aij = —c^{x^) ^ the inverse of A*^ exists. 
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7.2 The heat conduction equation with a flux in a Riemannian space 

The heat equation with a flux in an rt— dimensional Riemannian space with metric g.ij is: 

H {u) ^ q{t,x\u) (79) 

where 

H {u) := g'^^Uij - Vui - ut. 

The term q indicates that the system exchanges energy with the environment. In this case the Lie point 
symmetry vector is: 

X = r {X\u)d,+T^{x\u)du 

where a = t,i. For this equation we have: 

A"^0, A" = 0, A'^ = g'\B' = r{t,x'),B' ^lj{x,u) =q{t,x'',u) . 
For this PDF the symmetry conditions ([57l) - ([62|) become : 

ri = a{t,x')u + b{t,x') (80) 

e=e{t) (81) 

g^^{aijU + bij) - {a,iU + bj)T' - {a^tu + b^t) + Ag = ^*g,t + £_^q,k + m,u (82) 

- 'i9'''a^ + ar'^ - ^'.r + fF^, - AF'^ - (83) 

^es.fi'y- = (a - A)5jj. (84) 

The solution of the symmetry conditions is summarized in Theorem [9l The proof of the theorem is given in 
Appendix B. 

Theorem 9 The Lie point symmetries of the heat equation with flux i.e. 

g^-'uij — T^Ui — Uf = q (<, X, u) (85) 

in a n-dimensional Riemannian space with metric gij are constructed form the homothetic algebra of the metric 
as follows: 

a. y is a nongradient HV/KV. 
The Lie point symmetry is 

X= {2c2ijt + ci)dt+C2Y'd^ + {a{t)u + b{t,x))du (86) 
where a{t), b (t, x'') , q (t, x^ , uj must satisfy the constraint equation 

-atu + H (6) - {au + b) + aq- {2ij}C2qt + Clq)^ - C2g,jF' = 0. (87) 

b. ^ S'' is a gradient HV/KV. 
The Lie point symmetry is 

X=(^2i: J Tdt + ci^ dt + TS^'d^ + ( ^TtS + F (t)^ u + b{t, x)J 9„ 
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where F{t),T{t),b {t,x^) ,q{t, x*' , u) must satisfy the constraint equation 

- {^-^Tt^ + ^TttS - Ft^ u + H{b) + 

^TtS + F^u + b^q^.^+(^-~^TtS + F^q-(^2HjTdt + c,q^ - Tq^,S'\ (89) 
We apply theorem [S] for special forms of the function q (i, x, u). 

7.2.1 The homogeneous heat equation i.e. q{t,x,u) = 
In this case we have the following results 

Theorem 10 The Lie point symmetries of the homogeneous heat conduction equation in an dimensional 
Riemannian space 

g'^u,j - ru, -ut^Q (90) 

are constructed from the homothetic algebra of the metric gij as follows: 

(a) IfY^ is a nongradient HV/KV of the metric gij, the Lie point symmetry is 

X = {2'4,cit + ca) dt + ciY'd, + {a^u + b (t, x')) du (91) 

where ci,C2, ,ao are constants and b (t, a;') is a solution of the homogeneous heat equation, 
(c) LfY^ — S''' is a gradient HV/KV of the metric gij the Lie point symmetry is 

X = [csijt^ + C4< + C5)dt + {est + c4)S'd, + (-y ^ - y W^i + cg ) u9„ + 6 {t, x') du (92) 

where C3, C4, C5 are constants and b (t, a;*) is a solution of the homogeneous heat equation. 

In order to compare the above result with the existing results in the literature we consider the heat equation 
in a Euclidian space of dimension n. Then in Cartesian coordinates gij — 6ij, and the heat equation is: 

S'^Uij -ut=0. (93) 

The homothetic algebra of the space consists of the n gradient KVs di with generating functions x^ , the 
nongradient KVs Xjj which are the rotations and a gradient HV with gradient function H = ROr. According 
to theorem [lU] the Lie symmetries of the heat equation in the Euclidian n dimensional space are (we may take 

V-^i) 

X = [c3Vi' + (c4 + 2V'Ci)t + C5 + C2]at+[cir' + (c3t + C4)5*]a, + (94) 
+ \ [an + y S* + -^ni't - cg j u + 6 [t, a;')J (9„. 

This result agrees with the results of |3j p. 158. 

Next we consider the de Sitter spacetime (a four dimensional space of constant curvature and Lorentzian 
character) whose metric is: 

ds' = 



2 (-dr^ + dx'^ + dy'^ + dz^) 



(1 + f (-t2 + a;2 + j/2 + ;22)y 
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It is known that the homothetic algebra of this space consists of the ten KVs 

^- + ji]dy + [y-A 

^ ~ + if ) 

^ — ^ " X ^ ^^"^^ ^ ^^^^ ^ ^^^^ 

Xs = a;9r + t9:e , = + rdy , X7 = + rd^ , = yd^ - xdy 
Xg = zdj; - xdz , Xio = zdy - 2/9^ 

all of which are nongradient. According to Theorem [TUl the Lie symmetries of the heat equation in de Sitter 
space are 

10 

9t + ^ caXa + {aou + b (x, 

A=l 

From Theorem (TU] we have the following additional results. 

Corollary 11 The one dimensional homogenous heat equation admits a maximum number of seven Lie sym- 
metries (modulo a solution of the heat equation). 

Proof. The homothetic group of a 1-dimensional metric ds^ — (x) dx^ consists of one gradient KV (the 
■^^dx) and one gradient HV j g{x)dx dx)- According to theorem [TOl from the KV we have two Lie 
symmetries and from the gradient HV another two Lie symmetries. To these we have to add the two Lie 
symmetries X — OQudu + b (t, a;*) du and the trivial Lie symmetry dt where b (t, a;*) is a solution of the heat 
equation. ■ 

Corollary 12 The homogeneous heat equation in a space of constant curvature of dimension n has at most 
(n + 3) + {n — 1) (modulo a solution of the heat equation). 

Proof. A space of constant curvature of dimension n admits n + {n — 1) nongradient KVs To these we have 
to add the Lie symmetries X = cdt + a^udu + b (t, a;*) (?„. ■ 

Corollary 13 The heat conduction equation in a space of dimension n admits at most (n + 3) + 5 Lie 
symmetries ( modulo a solution of the heat equation ) and if this is the case the space is flat. 

Proof. The space with the maximum homothetic algebra is the flat space which admits n gradient KVs, 
in {n — 1) nongradient KVs and one gradient HV. Therefore from Case 1, of the theorem we have (n + 1) + 
in (n — 1) Lie symmetries. From Case 2. we have another (n + 1) Lie symmetries and to these we have to add 
the Lie symmetries X = cidt + OQudu + b (t, a;*) du where b (t, x*) is a solution of the heat equation. The set of 
all these symmetries is 1 + 2n + ^n (n - 1) + 2 + 1 + 1 = ^n (n + 3) + 5 [11]. ■ 
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7.2.2 Case q{t,x,u) ^ q{u) 

In this case we have the following result: 

Theorem 14 The Lie symmetries of the heat equation with conduction q(u) in a n dimensional Riemannian 
space 

g'^u^j -Tui-ut ^ q (u) (95) 

are constructed form the homothetic algebra of the metric as follows. 

a. is a HV/KV 
The Lie point symmetry is 



X = {2ctljt + ci)dt + cY'di + {a{t)u + b {t, x)) 

where the functions a (i) , b (t, x) and q (u) satisfy the condition 

— atu + H {b) — {an + b) q^u + (a — 2-00) q — Q. 

b. = S'' is a gradient HV/KV 
The Lie point symmetry is 

X^(^i, j Tdt + ci^dt+TS-'d,+ (^(^-^TtS + F{t)^u + bit,x)^du 



(96) 



(97) 



(98) 



where b {t, x) is a solution of the homogeneous heat equation, the functions T{t), F (t) and the flux q (u) satisfie 
the equation: 



lTt,p + ^TttS-Ft]u + Hib) 



-^TtS + F^u + b^ q,u + (-ItiS + F]q- 2i^qT = (99) 



For various cases of q (u) we obtain the results of the following tabic. 



q{u) 


Lie Symmetry vector 


qou 


{iPTot^ + 2ajjt + ci) dt + {cY' + TotS-') c 


%+ 




+ ( [-2yjcqot + aQ + To{-^S - ^Pqot^ - 


u + b{t, x)) du where H{b) - bqo ^ 




{2cipt + ci) dt + cY'd., + (j^u) du 




ulnu 


cidt + {Y' + Tae-^K') d, + {aoe-^u) du 




e" 


{2cil}t + ci)dt+ cY'd^ + {-2iPc) du 





where is a HV/KV , S*'' is a gradient HV/KV and K'^ is a gradient KV. 

8 Conclusion 

The main result of this work is Proposition [5] which states that the Lie symmetries of the PDEs of the form (|36) 
are obtained from the conformal Killing vectors of the metric defined by the coefficients Ay , provided Aij^u — 



17 



This result is quite general and covers many well known and important PDEs of Physics. The geometrization 
of the Lie symmetries and their association with the collineations of the metric, dissociates their determination 
from the dimension of the space, because the collineations of the metric depend (in general) on the type of the 
metric and not on the dimensions of the space where the metric resides. Furthermore, this association provides 
a wealth of results of Differential Geometry on collineations which is possible to be used in the determination 
of the Lie symmetries. 

We have applied the above theoretical results to two cases. The first case concerns the two dimensional 
wave equation in an inhomogeneous medium and shows the application of the general results in practice. The 
second example concerns the determination of the Lie symmetry vectors of the heat conduction equation with a 
flux in a Riemannian space, a problem which has not been considered before in the literature. We proved that 
the Lie symmetry algebra of this PDE is generated from the homothetic algebra of the metric. We specialized 
the equation to the homogeneous heat conduction equation and regained the existing results for the Newtonian 
case. Finally it can be shown that the Lie symmetries of the Poisson equation in a Riemannian space which 
have been computed in [12^ are obtained from the present formalism in a straightforward manner. 



Appendix A 



We prove the statement for 7i ~ 2. The generalization to any n is straightforward. For a general proof see |20) . 
We consider A^^ as a matrix and assume that the inverse of this matrix exists. We denote the inverse matrix 
with Bij and we get form ((4T|) : 







^'n = 0. (100) 

Now assume that the matrix A'^ does not have an inverse. Then we consider n = 2 and write: 

All Ai2 



[A'- 



^detA'^ =A''A'^ -{A' 







Al2 A22 

where at least one of the A''^ ^ 0. Assume A^^ ^ 0. Then equation (|41]) for i — j — k — 1 gives 

3A"e'„ = ^ e,u = 0. 

The same equation for i = j = k = 2 gives 

therefore either = or = 0. If A^^ ~ then from the condition det A'^ = we have A^"^ = hence 
Aij = which we do not assume. Therefore 0. 
We consider now equations i = j ^ k and find: 

Because i k this gives A"C''ti = and because we have assumed A^^ 7^ it follows Therefore again we 
find = 0- 
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9 Appendix B 

Condition (|84l) means that is a CKV of the metric gij with conformal factor a{t,x^) — X{t,x^). 

Condition 1^ implies ^'^ = T [t) (x^) where is a HV with conformal factor -0, that is, we have: 



and 



from which follow 



LyiQij — '^tpgij , i/'=constant. 



Condition (1821) becomes: 



(103) 



^* (t) = 2?/' y Tdi. (101) 
- 2g'''a^, + TtY'' = 0. (102) 



H{a)u + H{b) + {a- e*t)g = C*qj + T{t)Y\,k + m.u ^ 

H{a)u + H{b) - (au + b) + aq - {^'q)^t - T{t)Y\^k - 

H{a)u + H (6) - {au + b) q,,, + aq - (^2^/jq J Tdtj - Tq,,Y' = 0. 

We consider the following cases: 
Case 1 y*^ is a HV/KV 

From (fT02)) we have that Tt = O^T (t) = C2 and a,^ = ^ a (i, x'') = a{t). Then (fT03)) becomes 

-atU + H (5) - {au + b) q.^ + aq - {2'4>C2qt + ciq)^ - C2q.,Y' = (104) 

Case 2 is a gradient HV/KV, that is Y^ = S*'*^ 
From (|102l) wc have 

a(t,x'=) =-iTt5 + F(t). (105) 
Replacing in ()103p we find the constraint equation: 

= (^-^r^^ + ]^TuS ~ F?j u + H{b) + 

^TtS + F^u + b^q,,+ (^-^TtS + F^q-(^2yjqjTdt + ciq^ -Tq^,S'\ (106) 
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